Dissolution of traffic jam via additional local interactions 
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We use a cellular automata approach to numerically investigate traffic flow patterns on a single 
lane. The free-flow phase (F), the synchronized phase (S), and the jam phase (J) are observed and 
the transitions among them are studied as the vehicular density p is slowly varied. If p is decreased 
from well inside the J phase, the flux $ follows the lower branch of the hysteresis loop, implying 
that the adiabatic decrease of p is not an efficient way to put the system back into S or F phases. 
We propose a simple way to help the system to escape out of J phase, which is based on the local 
information of the velocities of downstream vehicles. 
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Since 1990s, traffic phenomena have drawn much at- 
tention amongst physicists 1]. Pure theoretical interest 
of the study of the traffic problem originates from its va- 
riety of fascinating dynamic emergent behaviors , like 
the phase transitions between different dynamical states, 
the existence of hysteresis behavior, and the historic de- 
pendence of traffic flow, to name a few. Recently, there 
has been a growing consensus on the existences of three 
different dynamic phases: when the vehicular density p is 
small enough, the system first stays in the free-flow (F) 
phase in which vehicles can have the maximum speed al- 
lowed, and thus the flux $ is linearly proportional to p. 
As p is increased, the speeds of vehicles are synchronized 
(S phase) to each other, without having stop-and-go mo- 
tion, and then eventually the jam (J) phase arises beyond 
some value of p. 

The cellular automata (CA) model proposed by Kerner 
et al. (KKW) [s*] to explain the three traffic phases, has 
been recently applied for multilane roads and compared 
with single vehicle data measured on American high- 
way [5| . Knospe et. al. Q extended Nagel-Schreckenberg 
model [6j| to include the anticipation effect and the re- 
duced acceleration capabilities, and a further extension 
was made by Jiang and Wu [3] with the difference of 
sensitivities between stopped and moving cars taken into 
account. Davis @ used the modifled optimal velocity 
model and also produced the synchronized flow phase, 
but with unrealistic speed of upstream front. In Ref. 0, 
realistic assumptions of the limited deceleration capabil- 
ity and human overreaction has been used as important 
ingredients of the CA model. Similar ideas of the limited 
decelerating capability and cautious driving when con- 
gested have been adopted in Ref. [ifl] , in which a further 
improvement was made by letting the maximum speed 
be dependent on headway distance to incorporate the 
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FIG. 1: (Color online) Fundamental diagram: the flux $ ver- 
sus the vehicular density p shows a hysteresis behavior. As 
p is increased from the free-flow (F) phase (a), <1> follows the 
upper branch of the curve, passing the synchronized (S) phase 
(b and c), and the jam (J) phase (d). When started from the 
J phase (d), $ follows the lower branch of the curve, passing 
the points e and f (both are in J phase, see Fig. [2]). The 
dotted blue curve in the middle is obtained when our method 
of dissolution of the jam by using an additional interaction is 
turned on. 



idea of the comfortable driving. A CA model with the 
similar deceleration limitation but without the slow-to- 
start effect has been proposed, which again reproduced 
the three traffic phases and discontinuous transitions be- 
tween them [Tlj , which in turn suggests that the slow-to- 
start effect is not a cause but a result of the transition 
to the wide-moving jam. In Ref. an attempt to in- 
terpret the traffic flow in the viewpoint of the complex 
network of traffic states was made. Very recently, the 
emergence of a jam without bottlenecks has been ob- 
served experimentally [l3l|, which clearly indicates that 
the traffic jam is a collective phenomenon resulting from 
interactions between vehicles. In the present paper, we 
follow the same line of reasoning and propose that one 
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can use an additional local interaction to dissolve the 
traffic jam. 

The main purpose of the present work is not to pro- 
pose another traffic model, but to study how the traffic 
jam can be dissolved through the use of the local infor- 
mation of car traffic. In this regard, we use the same CA 
model as in Ref. [9] , which has the mechanical restriction 
(implemented as limited acceleration and deceleration ca- 
pabilities) and the human overreaction as two important 
ingredients of the model. Throughout the present study, 
discretized one unit of time (distance) amounts to At =1 
s {Ax =1.5 m), and x^(i;Jj) denotes the position (the ve- 
locity) of the nth vehicle at time t. The human driving 
behavior is described by the two-state variable 7^ , which 
is set to 7* = for < vl^_^_i < t'^_|_2- We also assign 
7^ = for u^^_2 > Wmax — 1 with the speed limit Wmax- If 
neither condition above is satisfied, we set — 1 instead. 
In words, if the speed is successively nondecreasing in the 
downstream direction the driver expects that the traffic 
condition is improving. She can also be optimistic when 
the second leading vehicle is fast enough. The variable 
7 is not quenched but can change in time, and describes 
the driving behavior of each driver. It is to be noted that 
both 7 = and 7 = 1 represent human overreaction and 
only the type of overreaction is different (optimistic and 
defensive). 

We sketch here only the basic ideas of the model in 
Ref. a driver in the optimistic (defensive) state with 
7 = 0(1) leaves less (more) distance to the leading ve- 
hicle than that required for her safety. In comparison, 
an optimistic driver tries to avoid a collision only up to 
some time steps. By introducing optimistic and defen- 
sive behaviors, human factor is modeled as an overreac- 
tion to the local traffic condition. It is also assumed that 
vehicles can neither exceed the speed limit Wmax = 20 
(corresponding to 108 km/h), nor drive in the opposite 
direction (w < 0). As another important ingredient of the 
model, each vehicle can decrease their velocity stochasti- 
cally, to implement the random driving behavior in real- 
ity. All the parameter values are set identical to those in 
Ref. to mimic empirical observations. In Ref. [l3|, the 
critical vehicular density separating the free-flow phase 
the congested flow phase for the experimental setup of 
circularly moving cars was found to be about 25 vehi- 
cles/km as observed in real highway traffic, and the speed 
of the downstream jam front was measured as 20km/h. 
In the present study, for example, the downstream speed 
of the wide moving jam is found to be 18.4km/h, similar 
to the empirical measurement. 

We first report our simulation result [ij] for the fun- 
damental diagram. We use the single lane of the size 
L =40 000 (corresponding to 60 km) under the periodic 
boundary condition. As an initial condition we uniformly 
distribute N = 600 vehicles [p = 10 veh/km] and run the 
simulation by using the model in Ref. [9f . In order to ne- 
glect transient behaviors, we ignore the first T time steps 

and measure the average speed Wavg = (1/2^) St=T+i ^('•) 
for another T steps, where the spatial average v{t) = 



(l/^)Sn=i^ra- '^^^ traffic flux is then computed by 
$ = pva-vg- We confirm that T = 30 000 is long enough 
to ensure that the system approaches the steady state. 
The slow increase and the decrease of p are implemented 
as follows: in order to disturb the flow as little as possi- 
ble, the addition of a single vehicle is done at the position 
where the distance between two cars (x„+i — a;„) is max- 
imum. When a car is removed from the road, it is done 
for the slowest car, in order to make the jam dissolve 
more easily. 

Figure [1] exhibits the (global) fundamental diagram, 
the flux $ versus the vehicular density p. The upper 
branch of the curve is obtained when p is slowly in- 
creased. When p = 50(veh/km) is approached from be- 
low, we begin to decrease p as described above. The dot- 
ted blue curve marked as "Mechanism On" in Fig. [l]will 
be explained below. As is well-known, the fundamental 
diagram exhibits hysteresis behavior: $ follows different 
curves as p is increased and decreased. The existence of 
the hysteresis loop indicates that once the traffic jam is 
formed, it is quite difficult to dissolve the jam unless the 
vehicular density is significantly decreased. On the other 
hand, if we decrease p when the system is at point c in 
Fig. [1] we confirm that the system can easily go back 
to the F phase. In other words, in the present model, 
there is no hysteresis behavior between the F phase and 
the S phase, which in turn indicates that the transition 
between S and F is not discontinuous in our study. In 
comparison, some of existing CA model studies have pro- 
duced a continuous F-S transition as in the present 
work while others supported the discontinuous transi- 
tion [H]. 

For better understanding of the characteristics of the 
three different phases, the spatio-temporal traffic pat- 
terns are displayed in Fig. [21 each panel corresponds 
to six different points (a-f) in the fundamental diagram 
(Fig. [ij. An advantage of the traffic model in Ref. Q 
is that we can easily see the connection between driving 
behavior and the traffic ffow. In the F phase [a in Fig. [1] 
and Fig. [2ja)], almost all drivers become optimistic and 
drive at the maximum speed allowed. As p is increased, 
the inhomogeneous S phase (b) is clearly observed be- 
tween (a) the F and (c) the homogeneous S phases. A 
further increase of p eventually results in (d) the J phase. 
Interestingly, we observe that the downstream front of 
the jam and the downstream front of the synchronized 
region move at different speeds, as can be compared in 
Figs. [H(b) and (d) [17]. We emphasize that the traffic 
fiow patterns in Fig. [2] are obtained after long periods: 
under the periodic boundary condition adopted in the 
present study, once the jamming occurs, it is not sponta- 
neously dissolved either to the F phase [from (f) to (a)] 
or to the S phase [from (e) to (c)]. The measured speed 
(w 50 km/h) of the downstream front of the synchro- 
nized region in Fig. [UJb) is much higher than the corre- 
sponding speed (~ 18.4 km/h) for the jammed region in 
Fig- Hie) and (f). In comparison, the improved Kerner's 
three-phase model in Ref. [l^ has produced a realistic 
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FIG. 2: (Color online) Spatio-temporal {x,t) patterns of traffic flows at points a - f in Fig. [T] Dots with bright and the dark 
colors (green and red online, respectively) represent vehicles at optimistic and defensive states. As the vehicular density p is 
increased (see Fig. transitions occur from (a) the F phase to (c) the homogeneous S phase and then to (d) the J phase. 
Between the F and the homogeneous S phase, the inhomogeneous S phase in (b) is observed. As p is slowly decreased, the 
system does not escape out of the J phase [compare (e) and (c), and (f) and (a), respectively]. 



speed (5 cells per time unit corresponding to 8.5kni/h) 
of the downstream front of the synchronized How, while 
in Fig. 3 of Ref. (llj the corresponding speed is measured 
to be roughly about 32km/h. We remark that the some- 
how unrealistically high speed of the downstream front 
of the synchronized flow in Fig. [Ifb) could be an arti- 
fact of the present model. Likewise, the absence of the 



hysteresis behavior around S-F transition could also be 
model dependent. 

The main objective of this study is to suggest a simple 
way to dissolve the traffic jam within the limitation of 
the present CA model. If strong enforcement by the cen- 
tral authority is allowed, the most efficient way to have 
maximum traffic flow is to directly control the vehicular 
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FIG. 3: (Color online) The average speed v{t) versus the time t. We start from the J phase (at f in Figs. 1 and 2), and turn on 
our jam dissolving mechanism (see text) a.t t = 5000. (a)-(c) The first method based on the speed of the m-th leading vehicle. 
(d)-(f) The second method based on the speed of the vehicle at the geographic distance d. 



density p so that traffic flow is tuned toward the point 
a in Fig. [TJ In this iron-fist solution, more cars are al- 
lowed to enter the road ii p < pc ~ 20(veh/km), while 
all the entrances are blocked if p > pc- Although not 
realistic, this centralized way of maximizing traffic flow 
could be a very efficient way thanks to the absence of 
the hysteresis behavior between the F and the S phases. 
In other words, although the central authority fails to 
detect the maximum-flow condition p = pc and thus the 
synchronized flow already starts, a simple reduction of 
p will put the traffic back to the free-flow state. In the 
real-world traffic, however, it is not plausible that drivers 
are willing to accept the centralized control. We believe 
that the above consideration raises an interesting issue 
in the game theoretic situation: near p = p^, the flow is 
maximum and every driver will complain much if she is 
not allowed to enter the road which her leading vehicle 
just entered in a few seconds earlier. For her own sake, 
she is better off entering the road against the will of the 
central authority. However, in the public perspective, the 
overall traffic flow becomes worse by allowing her to get 
in. A similar situation has been termed as the price of 
anarchy in the study of urban car traffic [l^ . 

It is to be noted that in real traffic situations, feedback 
control of traffic flux at ramps is better suited to dissolve 
the traffic jam than the control of the global average ve- 
hicular density. The spatial inhomogeneity and the un- 



avoidable time lag between the start of the control and 
the system's response require adaptive control on ramps. 
In this regard, " ramp- metering" has already been well 
studied and applied in real highway systems [20l]. 

We suggest below a noncentralized way to dissolve the 
traffic jam which is based on local traffic information. 
Traffic flow can be of course improved by building a new 
road. However, the construction of new roads is not a 
practical alternative anymore in most countries pTj . Ac- 
cordingly, a recent focus towards better traffic flow has 
been made regarding how we can use existing roads in a 
more intelligent way. One of the promising directions for 
achieving better traffic flow is by using the so-called adap- 
tive cruise control [111, The basic idea behind this 
is that each car can be equipped with a detector which 
senses local traffic information and sends it to nearby 
cars. This gathered information can then be used to 
decide a driving strategy. In the same spirit, we pro- 
pose below a jam dissolving mechanism which is based 
on downstream traffic information. 

The robustness of the jam phase indicates that the 
incoming flux to the jam is not smaller than the outgo- 
ing flux. Accordingly, if an appropriate softening of the 
myopic behavior of jam-approaching drivers is made, it 
could reduce the incoming flux, resulting in the eventual 
dissolution of the jam. From this reasoning, we propose 
a simple way to dissolve the jam: if a vehicle is approach- 
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ing the upstream jam front, the driving behavior becomes 
defensive. In this regard, the use of the model in Ref. [3] 
has a great benefit since the model already contains the 
dynamic variable 7 which describes the driving behav- 
iors. Accordingly, our method is written as 

• If the mth leading vehicle has speed lower than the 
threshold value Vth, (i-e., if wJi+m < ''^th), the driver 
becomes defensive (7^ = !)• 

In the broad range of m and vth we indeed observe that 
the traffic jam is dissolved. We then try to optimize the 
values of m and Vth to make the dissolution of the jam 
happen as quickly as possible. As a result, we find that 
when Vth — 8(w 43 km/h) and 10 < m < 30 the jam is 
dissolved faster than the use of other values of vth and m. 
In Fig. |3Ia)-(c), we display the jam dissolving behavior 
for various values of m: even at to = 3 the jam is found 
to dissolve eventually, as shown in Fig. |3Jc). It is note- 
worthy that the congested traffic state is characterized by 
V <40 km/h in Ref. 2l|, similar to our optimized value 
of Vth ■ We also use a modified version of our method in 
which the geographic distance d to the downstream pre- 
ceding vehicle is used instead of m: If the farthest vehicle 
in the downstream direction at a distance shorter than d 
has a speed lower than vth , 7 is set to one (defensive driv- 
ing). We find that the jam is dissolved in a broad range 
of d for c? > 50 corresponding to d > 75 m [see Fig. [S^d)- 
(f )] . We also observe that if the fraction of vehicles with 
our method turned on is larger than about 80%, the jam 
is eventually dissolved for both types of our methods. 

Figure [4] exhibits the dissolution of the jam in the 
spatio-temporal plot {x-t) [(a) and (b)], and (c) the 
spatial average v{t) as a function of time t, at p = 
20.0(veh/km). The initial state at t = is obtained from 
the long-time simulation at point f in Figs. [T] and [H and 
the mechanism explained above is on at t —5 000 with the 
parameters m = 15 and Vth — S- It is clearly seen that as 
time goes on, more and more jam-approaching drivers 
become defensive [Fig. [H^a)], and the jammed region 
with higher local density becomes narrower [Fig.UJ^b)], as 
the upstream jam front eventually catches up the down- 
stream jam front at t ^ 14000 [Fig. [4ja)]. Interestingly, 
a narrow track of defensive drivers persists for a quite 
long time, until all drivers become optimistic at around 
t = 18 000. In Fig. [1] we also include the curve when 
our method of using additional downstream traffic infor- 
mation is turned on at all times. The key observations 
in Fig. [T] are as follows: (i) The hysteresis behavior dis- 
appears and thus increase and decrease of p result in the 
same curve in the fundamental diagram, (ii) The use of 



our method greatly improves the traffic flow in a broad 
range of vehicular densities for p < 26, and effectively 
removes the jam phase in the lower branch, (iii) When 
26 ^ p ^ 31, our method still helps but the traffic flow 
is worse than the upper branch without it. (iv) Inter- 
estingly, when p 31, the traffic flow with our method 
on is worse than the jam phase without it. (v) However, 
the jam that is robust in the broad range of 20 < p < 30 
is efficiently dissolved by applying our simple method, 
which leads to a considerable flux increase compared to 
that of the jam phase in this density region. 

For comparisons, we also display in Fig. [5] the spatio- 
temporal patterns and the average velocities versus time 
at p = 30 [(a) and (b)], and p = 45 [(c) and (d)], re- 
spectively. The mechanism is turned on at f = 5 000 
as in Fig. 01 At p = 30, the congested region eventu- 
ally spreads out and the average velocity in the steady 
state becomes larger than the value at the lower branch 
in Fig. [T] [see Fig. HKb)]. In contrast, when the vehicular 
density is too high as in Fig.[5l^c), the steady state value 
of the average velocity is smaller than the one with the 
mechanism off [see Fig. EKd) and (iv) above] . 

In summary, we have used the single-lane traffic model 
proposed in Ref. ,9i| to study the nature of three traffic 
phases. As the vehicular density is slowly increased and 
decreased, hysteresis behavior is observed in the funda- 
mental diagram for the flux and the density. Within the 
limitations of the used CA model, we propose a simple 
method to dissolve the jam, in which driving behavior 
is forced to become defensive if the leading vehicle a few 
cars ahead is dangerously slow. We believe that the main 
conclusion of jam dissolution via slowing down the jam- 
approaching vehicles by informing them of the speeds of 
leading vehicles could be general. However, the valid- 
ity of our results needs to be checked more carefully in 
further studies. Our proposed mechanism dissolves the 
jam since it reduces the inflow without altering the out- 
flow from the jam. Combined with the complementary 
method of increasing the outflow by using automated 
cruise control, the elimination of the jam can be made 
more effective. 
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dots) in the intermediate times, and eventually the traffic jam is dissolved for all drivers become optimistic (green dots), (b) 
The expansion of a part of (a) to show more clearly that the jammed region becomes narrower after t — 5000. (c) The spatial 
average speed v{t) versus time t. 
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p = 30 veh/km p = 45 veh/km 




FIG. 5: (Color online) Spatio-temporal {x,t) pattern of the traffic flow starting from the J phase at (a) p — 30 and (c) p = 45. 
The average velocity v{t) as a function of time t for the corresponding spatio-temporal pattern is also shown in (b) and (d). It 
is to be noted that at p = 30, our jam dissolving mechanism improves the flow, while at p = 45 it makes the flow worse (see 
Fig.0. 



